Axial anomaly as a collective effect of meson spectrum 
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Abstract 

We study the transition form factors of the light mesons in the kinematics, where one photon 
is real and other is virtual. Using the dispersive approach to axial anomaly we show that the 
axial anomaly in this case reveals itself as a collective effect of meson spectrum. This allows us 
to get the relation between possible corrections to continuum and to lower states within QCD 
method which does not rely on factorization hypothesis. We show, relying on the recent data 
of the BaBar Collaboration, that the relative correction to continuum is quite small, and small 
correction to continuum can dramatically change the pion form factor. 

1 Introduction 

The phenomenon of axial anomaly [H[2] is known to be one of the most subtle effects of quantum 
field theory. Perhaps the most vivid manifestation of it in particle physics can be found in two- 
photon decays of pseudoscalars. Usually the differential form of the axial anomaly is utilized to 
study this kind of processes. However, less known dispersive approach to axial anomaly ( [3j|4], f° r 
a review, see [5]) leads to anomaly sum rule (ASR) relation, which provides a very powerful tool for 
study of various characteristics of meson spectrum. The absence of corrections to ASR allows to 
get certain exact relations between characteristics of hadrons, such as decay width [6] and relations 
between mixing parameters of pseudoscalars [7JE]. 

In the study of two-photon decays of pseudoscalars, usually the case of real photons is considered. 
However, ASR is valid for virtual photons also [9 |I10| which leads to several interesting applications. 
As we will see, in the case of one real and one virtual photon [9] the ASR takes into account the 
infinite number of meson states (which can contribute to it due to their quantum numbers), i.e. 
axial anomaly reveals itself as a collective effect of meson spectrum. This nontrivial fact together 
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with the exactness of ASR gives us a tool to study different characteristics of meson spectrum (form 
factors of mesons, relation between possible corrections to lower meson states and continuum). 

The experimental measurements of the photon-pion 77* — > ir° transition form factor F n ^(Q 2 ) 
in the range of photon virtualities Q 2 < 9 GeV 2 were performed by CELLO [11] and CLEO |12j 
Collaborations. The measured values of F 7rry (Q 2 ) appeared to be consistent with predictions based 
on factorization approach to pQCD. Surprisingly, recent data of BaBar Collaboration on F^{Q 2 ) 
|13j . which is available in the range 4 < Q 2 < 40 GeV 2 , have shown a strong disagreement with 
the pQCD predicted behaviour of 77* — > it transition form factor. 

Though the BaBar data in the range Q 2 < 10 GeV 2 fit well the CLEO data and is in a good 
agreement with theoretical predictions from the light-cone QCD sum rules (LCSR), offered in [14], 
however, at larger virtualities strong disagreement takes place. Moreover, more precise recent 
LCSR analysis |15U16j shows, that it is impossible to explain BaBar data on F n j(Q 2 ) at large Q 2 
by use of usual (endpoint-suppressed) form of pion distribution amplitude. This have led (quite 
unexpectedly) to the question of pQCD factorization validity. Recently, there were proposed several 
approaches to explain such anomalous behaviour of F nl (Q 2 ) [T7H2T] . in particular, questioning 
pQCD factorization. At the same time, in [22] the authors give some arguments against the related 
flat-type pion distribution amplitude, while in [23] some doubts about BaBar results analysis were 
expressed. 

In this paper we study what can be learnt about the meson-photon transition form factors from 
the anomaly sum rule for the case of one virtual photon. This generalizes the usual application of 
anomaly, providing the boundary condition in the real photon limit only. Our (non-perturbative) 
QCD method does not imply the QCD factorization and is valid even if the QCD factorization 
is broken. It is shown, that using axial anomaly in the dispersive approach we can get the exact 
relations between possible corrections to lower states and continuum providing a possibility of 
relatively large corrections to the lower states. 

2 Anomaly sum rule 

Following [9], we briefly remind some results of the dispersive approach to axial anomaly which are 
relevant for this paper. The VVA triangle graph amplitude 

T a ^(k,q) = J d 4 xd 4 ye^ + ^)(0|T{j5(0)J / ,(x)J,(y)}|0) (1) 

contains axial current = (uj^ a u — d"f^ a d) and two vector currents J p = ((2/3)«7 M n — 
{l/2>)d^ p d)\ k,q are momenta of photons. This amplitude can be presented as a tensor decom- 
position 

T a [n/(k, q) — F\ s apiVp k p + F2 SafivpQ^ 

+ F 3 q v e app(J k p q a + F A q v e app(J k p q a (2) 
+ F 5 k p e aupa k p q a + F e q p e aupa k p q u , 

where the coefficients Fj = Fj (k 2 , q 2 , p 2 ; m 2 ) , p = k + q, j = 1, . . . , 6 are the corresponding Lorentz 
invariant amplitudes (form factors). Note that these form factors do not have kinematical singu- 
larities and are suitable for dispersive approach, which we use to derive anomaly sum rule. 
Symmetries of the amplitude T apu (k,q) impose the relations for the form factors Fj(k,q). 
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Bose symmetry, i.e. T afll/ (k,q) = T au ^(q,k) leads to 

F 1 (k,q) = -F 2 (q,k), 

F 3 (k,q)=-F 6 (q,k), (3) 
F i (k,q) = -F 5 (q,k). 

Vector Ward identities for the amplitude T afJil/ (k,q) 

WT atlv = 0, q v T aflP = (4) 

in terms of form factors read: 

F 1 =k.qF 3 +q 2 F 4 , F 2 = k 2 F 5 + k.q F & . (5) 
Anomalous axial- vector Ward identity for T Qlll/ (k,q) [TJG] 

p a T afll/ {k, q) = 2mT^ v {k, q) + -L e ^ p(T Fg CT (6) 

in terms of form factors can be rewritten as follows: 

F 2 - F x = 2mG + » (7) 

where G is a form factor, related to the 2nd rank pseudotensor Tn V involved in the "normal term" 
on the r.h.s. of ([U]): 

V(M) = G£„^T. (8) 

Writing the unsubtracted dispersion relations for the form factors one gets the finite subtraction 
for axial current divergence resulting in the anomaly sum rule which for the kinematical configu- 
ration we are interested in (k 2 = 0, q 2 ^ 0) takes the form [9]: 

A 3a (t;q 2 ,m 2 )dt = — , (9) 
4m 2 2tt 

where 

A 3a = hm(F 3 -F 6 ). (10) 

It holds for an arbitrary quark mass m and for any q 2 in the considered region. Another important 
property of the above relation is absence of any a s corrections to the integral [23]. Moreover, it is 
expected that it does not have any nonperturbative corrections too ('t Hooft's principle). 

3 Transition form factors of mesons 

The form factor of the transition ir° — > 77* is defined from the matrix element: 

J d*xjtefrO(p)\ T {J^x)J v (p)}\0) = e^WtfF^ , (11) 

where k,q are momenta of virtual photons, p = k + q, and = ((2/3)^7^ — (1/3)^7^) is the 
electromagnetic current of light quarks. 
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Three-point correlation function T afll/ (k,q) has pion (pole at p 2 = m 2 ) and higher states con- 
tributions: 

T a u,u(k, q) = p a k p q a e llupcr F 7T1 + (higher states) , (12) 

p z — m z 

where f n is a pion decay constant, which can be defined as a coefficient in the projection of axial 
current J\ onto one-pion state: 

(0\4(0)\Tr°(p))=iV2p a U. (13) 

The pion decay constant /„- = 130.7 MeV is experimentally well determined from the decay of 
charged pion ir~ — > [i~v. Using the kinematical identities 

we can single out the pion contribution to ^(^3 — Fq). Then, the contribution of pion to Im(F 3 — 
Fq)/2 (imaginary part is taken w.r.t. p 2 ) is: 

^Im(F 3 - F 6 ) = V2f^F^(Q 2 )5{s - ml) , (15) 

where Q 2 = -q 2 0. 

It is well known that at Q 2 = the pion contribution saturates anomaly sum rule ([9]) [3] and 
F™ is known to be normalized by the 7r° — > 77 decay rate [2]: 

F - (0) = i^T (16) 

On the other hand at Q 2 7^ 0, factorization approach to perturbative quantum chromodynamics 
(pQCD) for exclusive process in the leading order in the strong coupling constant predicts [251(26] : 



FiryiQ ) = / ' + 0{l/Q^ , (17) 

where f n = 130.7 MeV and (p-x(x) is a pion distribution amplitude (DA). The pion DA depends 
on the renormalization scale [251127] and at large Q 2 asymptotically acquires a simple form [28j 
^asymp ^ = Qx (1 — x). This leads to asymptotic behaviour for the pion form factor: 

F^ mp (Q 2 ) = ^ + 0(1/Q 4 ) . (18) 
From (|17p and f)15|) we get the contribution of pion to anomaly sum rule (J9j) : 

2nf 2 /Q 2 . (19) 

We see, that at Q 2 7^ anomaly sum rule © cannot be saturated by pion contribution H due to 
1/Q 2 behavior, so we need to consider higher states. The higher mass pseudoscalar states have the 
same behavior and suppressed by the factor m 2 /m 2 es as follows from the PCAC (since d^J^ should 



1 Note that in the case of the choice of the basis differing from @ for AVV amplitude the results for transition 
formfactors do not change. 

2 This may be compared with Q 2 dependence [29| of pion-to-photons matrix element emerging in the analysis of 
sum rule [291130] for photon spin structure function. 
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vanish in the chiral limit). The other contributions are provided by axial mesons, the lightest of 
which is the ai(1260) meson. In fact, the contribution of longitudinally polarized a\ is given by the 
similar equation to (|18p at large Q 2 . indeed, the analysis of axial current bilocal matrix elements 
of oi is completely similar for that |31| of vector current matrix elements of p. The contribution of 
transversally polarized at to ([9]) decreases even faster. Actually, the same is true for all the higher 
axial mesons and mesons with higher spin. So we make an important observation: for the case 
Q 2 7^ the anomaly relation ([9]) cannot be explained in terms of any finite number of mesons due 
to the fact that all transition form factors are decreasing functions. That is why we conclude that 
only infinite number of higher states can saturate anomaly sum rule and therefore at Q 2 7^ the 
axial anomaly is a genuine collective effect of meson spectrum in contrast with the case of two real 
photons Q 2 = 0, where the anomaly sum rule is saturated by pion contribution only. Let us note 
that this conclusion does not depend on any choice of meson distribution amplitudes. 



4 Quark-hadron duality 

Now we proceed to particular analysis of anomaly sum rule ©using quark-hadron duality. In some 
sense the discussion of the previous section is based on the quark-hadron duality. In this section 
we will apply local quark-hadron duality to anomaly sum rule ([9]). According to the quark-hadron 
duality, let us saturate the the spectral density A^ a by pion and continuum contributions: 

A 3a (s, Q 2 ) = V2tt U5(s - ml)F^ (Q 2 ) + A^ D 9(s - s ), (20) 

where continuum contribution A^' D 9(s — so) as usually supposed to be equal to QCD calculated 
spectral function A^ a , while so is a continuum threshold. 

Substituting (|20p into ([9]), we obtain the anomaly sum rule in the following form: 

— = v^tt/^CQ 2 ) + J dsA^ D . (21) 

One-loop PT calculation [32|,l33j leads to a simple result for the spectral density function: 

^ Ci W>^(ITW> (22) 

so we can rewrite the anomaly sum rule in the following form: 

1 1 r°° O 2 

S = ^^«"> + SJL *(TfW- (23) 

and finally the pion form factor is 

F - (Q2 >=vfeV (24) 

where sq = 0.7 GeV 2 is a continuum threshold. 

This result coincides with the interpolation formula proposed by S. Brodsky and G.P. Lepage [26] 
(«o = 4vr 2 /2) 

F ^ {Q2) = 2V2n 2 f 7r l + Q 2 /(i7r 2 f 2 ) ' (25) 
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Figure 1: (Colour online). Relative contributions of 7r° (thick blue curve), a\ (thin blue curve) 
mesons (intervals of duality are 0.7 GeV 2 and 1.8 GeV 2 respectively), and continuum ( dot-dashed 
black curve) (continuum threshold is S\ = 2.5 GeV 2 ) to the anomaly sum rule (dashed red line). 

which was derived in the quark-hadron duality context by A.V. Radyushkin |34| . The contact with 
exact anomaly sum rule observed here allows to find the relations between contributions of ir° with 
higher states which may be chosen as a% and continuum. With account of these contributions the 
anomaly sum rule © can be rewritten in the following form: 

1 1 f°° G 2 

- = V2*f,F,^) +/„ + - £ dsj^ , (26) 

where I ai is a contribution of a\ meson to sum rule(which can be expressed in terms of a\ form 
factors), s\ = 2.5 GeV 2 is a continuum threshold for this case. 

Using the asymptotic formula for pion form factor (|24j) we can estimate the behavior of I ai at 
large Q 2 as 

/ - 1 O 2 Sl ~ S ° (27) 
ai 2vr V ( Sl + Q 2 )(so + Q 2 ) ' { ] 

This equation can be treated as a good interpolation for a\ contribution with correct asymptotic 
behavior (large and small Q 2 ). 

The plot for contributions of pion, a\ meson and continuum is shown in Fig.l. The figure illus- 
trates the anomaly collective effect: indeed, the contribution of infinite number of higher resonances 
(continuum contribution) dominates starting from relatively small Q 2 ~ 1.5 GeV 2 . 

5 Corrections' interplay and experimental data 

As we learned above, axial anomaly is a collective effect of meson spectrum. That is why it is 
natural study the relation between possible meson and continuum corrections from anomaly sum 
rule. Let us write out anomaly sum rule once more in the following form: 
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1 f°° 

~ — = / A 3a ( S ; Q 2 )ds = h + I ai + leant , (28) 

2vr Jo 

where continuum contribution I con t takes into account all other higher mass axial and higher spin 
states. 

As we already mentioned, the anomaly sum rule ([9]) is an exact relation (J °° A 3a (s; Q 2 )ds does 
not have any corrections). However, the continuum contribution 

/oo 
A 3a (s;Q 2 )ds (29) 
-i 

may have perturbative as well as power corrections. 

Note that the two-loop corrections to the whole triangle graph were found to be zero [35] 
implying the zero corrections to all spectral densities^. To match this result with the earlier found 
non-zero corrections in the factorization approach (see [16] and references therein) one should be 
careful. When one is dealing with the corrections to the form factor itself, only the corrections to the 
coefficient function should be considered, while all other corrections are absorbed to the definition 
of the distribution amplitude. At the same time, when the factorization theorem is applied for 
calculations of the large Q 2 asymptotics of VVA diagram, all the corrections should be taken into 
account. To do so, one should add Eq. (3.11) of [16j with the projector to the local axial current 
(proportional to asymptotic pion distribution amplitud^) and Eq. (Bl) (coinciding with Eq. (1) 
of the Erratum) of [38] and get a zero result, compatible with |35j . 

Therefore, the model of the corrections to continuum discussed below should rather correspond 
to some non-perturbative corrections. Let us first consider the contributions of local condensates. 
Naively, they should strongly decrease with Q 2 compensating the mass dimension of gluon (as 
quark one is suppressed even more) condensate. However the 't Hooft's principle requires (see [9], 
Section 4) the rapid decrease of the corrections with Borel parameter M 2 (related to s) so that 
the power of Q 2 in the denominator may be not so large. In reality the actual calculations do 
not satisfy this property and the situation may be improved by the use of non-local condensates 
(see [9] and references therein). Another possibility is other non-perturbative contributions, like 
instanton-induced ones. So we assume the appearance of such corrections in what follows modelling 
the corrections to continuum. 

In order to preserve the sum rule ()28f) (or in particular (j23[) . (I26P ) the corrections to continuum 
contribution should be exactly compensated by corrections to lower states, in particular to pion. 
It turns out that this is a rather uncommon situation: corrections to continuum are compensated 
by the main terms of lower states which are of the same order in Q 2 as continuum corrections. 

To be more specific, let us consider the model "7r°-(-continuum". From (|23p the main contribu- 
tions of pion and continuum read: 

1° = V2irUF^(Q 2 ) = ^-—^ , (30) 

Z7T Sq + Q 

Icont = 7T~ . 7v) " (^) 
27T S + Q 



3 Zero two loop corrections to spectral densities were also found in the massive case |36) . although later the 
inconsistency of this result with asymptotic mass expansion was pointed out [37] . 
4 We are indebted to S.V. Mikhailov for clarification of this point. 
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If the corrections to pion and continuum to ASR are 81^ and 5I CO nt respectively 

I v = 7° + 5I„ , (32) 

Icont = Icont ^Icont (33) 

then, since dl n = —SI CO nt, the ratio of relative corrections to continuum and pion is 

5 Icont/ Icont 



For instance, for Q 2 = 20 GeV 2 , s = 0.7 GeV 2 the ratio is 

SIcont/ Icont 



0.03 . (35) 



We see, that the relative correction to continuum is suppressed by factor 1/Q 2 as compared to the 
correction to pion. To illustrate our conclusion, we assume the correction to continuum at large 
Q 2 is 5I C ont = — cso ln ^ g^ +b . This correction preserves asymptotics of continuum contribution at 
large Q 2 . Contributions of pion and continuum to ASR then have the following explicit forms: 

1 Q 2 ln(Q 2 /s ) + b 

cs , (36) 



l cont 



2ir s + Q 2 u Q 



I * = KZ „ i 02 + cg o 772 • ( 3? ) 



i *o HQ 2 /s o) + b 
2vr s + Q 2 CSo Q 2 

If this correction corresponds, as it was discussed above, to (non-local) gluon condensate (G 2 ) one 
may formally substitute the dimensional factor cso by c{G 2 ) / sq stemming from (G 2 )/M 2 for Borel 
transforms. 

One can see, that the leading power correction to continuum preserving its asymptotics results 
in a substantial (of the order of the main term I®) contribution to the pion state changing the pion 
form factor asymptotics at large Q 2 . 

The experimental data on pion form factor behavior at large Q 2 allows us to get estimation for 
the corrections to continuum. If the pion form factor expression ()24j) have matched the experimental 
data, the continuum leading correction might be only of order l/Q 4 . However, the last BaBar data 
on pion transition form factor [13] manifests large discrepancy of F n ^ values at large Q 2 with the 
expected asymptotic behaviour (j24j) . 

Relying on the BaBar data, we can fit parameters b, c: 

b = -2.74, c = 0.045. (38) 

The corresponding plot of combination Q 2 F 7T ~ f for the best-fit parameters (I38D is shown in Fig. 2. 

Basing on the model "pion+continuum" we can calculate the relative correction to continuum 
contribution to ASR 5I CO nt/ 'Icont relying on different fits of F*^: 

\SIcontZlU = ^UF;, - F^)/Il nt = 2 ^ f p + Q2) Kr ~ |r (39) 

In Fig. 3 the ratios 8I CO nt/ Icont f° r our fit l38l) and fit, obtained in recent paper [18] are shown. 
We see, that the correction to continuum is indeed small, even though the BaBar data shows that 
the relative correction to pion contribution is large. 
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Figure 2: Q 2 F 1T1 {Q 2 ): BaBar data |13| . pQCD calculated behaviour (j2l|) (dashed curve) and 
logarithmically enhanced behaviour due to small correction to continuum (|371 138 j) (solid curve). 
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Figure 3: (Colour online). Relative correction to continuum 5I con t/Ic 0nt (Q 2 ): our fit (|36l l38j) (solid 
curve) and fit, obtained in [18] (dashed blue curve). 
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Above we estimated the correction to continuum in the model "7r°+continuum". However, one 
can consider more refined models like "7r° + ai+continuum". This model may assume interplay 
between corrections to three terms in the ASR. 

Moreover, in the case of the small correction to continuum which currently seems to be the 
most likely situation, ASR will lead to the relation between the transition form factors of pion and 
ai which may be further studied both theoretically and, most important, experimentally. 



6 Discussion and Conclusions 

Dispersive approach to axial anomaly proved to be a useful tool for studying the properties of 
meson spectrum. It is well known, that when both photons are real (Q 2 = 0) the ASR saturates 
by pion contribution only. However, when one of the photons is virtual (Q 2 7^ 0) we immediately 
get different situation: ASR can be saturated only with a full meson spectrum (any finite number 
of mesons cannot saturate the anomaly sum rule). So the axial anomaly is a collective effect of 
meson spectrum. 

The anomaly sum rule and quark-hadron duality in case of model "7r°+continuum" allows 
to reproduce the well-known Brodsky-Lepage interpolation formula for F w ~. We estimate the 
contribution of a\ meson to anomaly sum rule in the model "7r°+ai+continuum". 

The exactness of the anomaly sum rule leads to the relation between corrections to continuum 
and lower mass states contributions. The last experimental data on pion transition form factor F n ^ 
at large Q 2 allows us to estimate the possible continuum correction. 

One can also consider 7*7 — > r] transition form factor in the same way. Considering the similar 
ratio for the relative contributions of rj meson and continuum as (I34p (with the continuum threshold 



one for 7r°: 



Sq ~ 2.5 GeV 2 ), we can estimate the relative correction to 77 to be several times smaller than the 



5I V /I° 



^ ~ 0.3. (40) 
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